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Abstract. As a canonical model for wetting far from thermal equilibrium we study a 
Kardar-Parisi-Zhang interface growing on top of a hard-core substrate. Depending on 
the average growth velocity the model exhibits a non-equilibrium wetting transition 
which is characterized by an additional surface critical exponent 8. Simulating the 
single-step model in one spatial dimension we provide accurate numerical estimates 
for 9 and investigate the distribution of contact points between the substrate and the 
interface as a function of time. Moreover, we study the influence of finite-size effects, in 
particular the time needed until a finite substrate is completely covered by the wetting 
layer for the first time. 



1. Introduction 

When a chemically inert surface is exposed to a gas phase particles may preferentially 
attach to the surface, forming a wetting layer of a different phase. By changing physical 
parameters such as pressure or temperature these systems may undergo a wetting 
transition at which the thickness of the layer diverges and becomes macroscopic. In 
many experimental situations it is reasonable to assume that the wetting layer is in 
thermal equilibrium with the substrate so that methods of equilibrium statistical physics 
can be applied [1,2]. More recently, however, wetting phenomena far from equilibrium 
attracted considerable attention [3-9]. These studies are motivated by the question 
whether wetting processes under non-equilibrium conditions, for example caused by a 
flux of particles or energy, may show physical properties that are qualitatively different 
from those of equilibrium wetting, especially at or in the vicinity of a wetting transition. 

In order to describe a wetting process far from equilibrium one has to resort to 
a dynamical description of the microscopic processes. This can be done either by 
introducing explicit solid-on-solid growth models [4,5] where the microscopical dynamics 
is described by appropriate rules for deposition and absorption of adatoms over the 
substrate, or by studying appropriate Langevin equations for non-equilibrium growth [8]. 
Concerning their universal properties all these studies were so far related to the Kardar- 
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Table 1. Critical exponents of wetting transitions in one spatial dimension 



Parisi-Zhang equation [10] in one spatial dimension with an additional potential term [3] 
^jM = VQ + aV 2 h ft t) _ dV d ^ ] + \(Vh(x, t)f + C (x, t) , (1) 

where h(x, t) is the height of the interface defining the thickness of the growing layer, 
vq, o and A are real parameters and ({x,t) is a white Gaussian noise. Here V(h) is 
a potential that takes the influence of the substrate into account. The nonequilibrium 
properties of this equation are caused by the non-linear term X(Vh) 2 which can be 
shown to be a relevant perturbation in the renormalization group sense that violates 
detailed balance. Obviously this term breaks reflection symmetry h — > —h so that in 
combination with a symmetry-breaking substrate the sign of A is expected to play a 
significant role. In the special case A = 0, however, it can be shown that the dynamics 
in a stationary state obeys detailed balance and is thus at thermal equilibrium. Hence 
by varying A we can study the crossover from equilibrium to nonequilibrium wetting. 

Without a substrate, i.e., for V = 0, a one- dimensional interface evolving according 
to the KPZ equation is known to roughen as w(t) ~ t 13 , where w(t) is the width (defined 
as the standard deviation of the heights) and (3 is the so-called growth exponent. In a 
finite system of size L this roughening process continues until the width saturates at 
w s ~ L a when the intrinsic correlation length £ ~ t 1//z reaches the system size. Here 
z and a = z{3 are the dynamical exponent and the roughness exponent, respectively. 
In the case of KPZ growth in one spatial dimension, where a fluctuation-dissipation 
theorem holds, the exponents are given by simple rational values [11], see Tabled 

In order to study nonequilibrium wetting, one has to introduce a substrate by 
imposing a boundary in the microscopic dynamics or by choosing an appropriate form 
of the potential appearing in Eq. (JTJ). The simplest way to introduce such a boundary 
is to impose the restriction h(x, t) > 0, corresponding to an infinite-step potential of the 
form 

V(h) = i° . (2) 

v ; \ oo if h< K J 

In order to avoid singularities in the continuous formulation, this infinite-step potential 
is usually replaced in Eq. ((TJ) by an exponential function 



V(h) = exp(-/i) 



(3) 
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Figure 1. Upper row: Typical spatio-temporal pattern of contact points 
(pinned sites) between interface and substrate in the SSM (see text) at 
criticality, starting with a flat interface at t = 0. Lower row: Final configuration 
of the interface at t = 5000. Note that for A < there are only few contact 
points, while for A > the density of contact points is much higher. 



that was found to display essentially the same physical properties. 

Models defined by the addition of this simple boundary to a KPZ roughening 
surface are known to exhibit a continuous wetting transition. In an infinite system 
the critical point of this transition is determined by the asymptotic average velocity 
t>oo = Hindoo t>(t) of a free interface. For > the interface moves away from the 
substrate so that it no longer influences the dynamics. For < 0, however, the 
interface is driven towards the substrate, approaching a steady state of finite roughness 
that is characterized by a certain height profile. At the transition point, where = 0, a 
scale-invariant behavior is observed. Obviously, the order parameter that characterizes 
this transition is the density of sites where the interface touches the substrate 



or - equivalently - the spatial average exp(— h(x, t)) in the continuous formulation. 

From a theoretical point of view the influence of a hard-core substrate without 
additional short-range forces can be interpreted as a boundary condition in the space 
spanned by x, t, and h. As in other critical phenomena, this boundary condition does 
not change the underlying universality class and the associated bulk exponents, meaning 
that the process still belongs to the KPZ universality class. However, the substrate as an 
additional boundary condition adds a new feature, characterized by an order parameter 
p s (t) which is associated with additional surface exponent 8. This order parameter was 
found decrease with increasing time as a power law [3] 
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The exponent 6 is expected to take three different values depending on the sign of 



Nonequilibrium wetting of finite samples 



4 



the parameter A of the nonlinear term. In fact, already a visual inspection of typical 
interfaces at criticality reveals significant differences depending on the sign of A. As can 
be seen in Fig. ^ for A < only few contact points remain, leading to a fast decay of 
p s (t), while for A > the density of contact points decays only slowly. 

Quite interestingly, it has recently been discovered, both through numerical 
simulations and analytical arguments [12,13], that also the complete synchronization 
phase transition (CST) in most spatially extended chaotic dynamical systems belongs 
to the same universality class of a bounded KPZ equation with a negative nonlinear term 
A. It can be shown that, when the local dynamic is smooth enough (more precisely: 
the linearization of the system around the synchronized state correctly describes the 
propagation of finite perturbations), the critical behavior of the CST is described by 
the so called multiplicative noise (MN) Langevin equation 



where the order parameter n(x, t) > is the coarse grained density of unsynchronized 
regions. Here a and D > are real parameters, and ({x,t) is t ne same noise as in 
Eq. (^). By a Cole-Hopf transformation [3] 



Eq. © turns exactly into the bounded KPZ equation ([TJ with A = —a = —D and 
an exponential potential as in Eq. (JHJ). This indicates that the MN equation is just a 
realization of the KPZ universality class and its non-trivial surface critical behavior. In 
particular the order parameter n(x, t) in Eq. © scales exactly in the same way as the 
density of contact points p s (t) and therefore its spatial average will decay as t~ e . 

Eq. (jnj) has been studied numerically in a series of papers [15-17] and its 
transition is said to belong to the multiplicative noise I (MN1) class. Moreover, the 
related multiplicative Langevin equation which can be derived via an inverse Cole-Hopf 
transformation from the bounded KPZ equation (0) with A > has been numerically 
analyzed in Ref. [18]. The associated transition is said to belong to the multiplicative 
noise II (MN2) class. 

Presently, no exact analytical derivation of the surface exponent 9 in the 
nonequilibrium case is known. While in Ref. [19] the surface exponent has been obtained 
for both the MN1 and MN2 classes with a mean-field like approximation suited for 
one spatial dimension, the presence of a strong-coupling fixed point and of essential 
divergences makes the exact values of the surface exponent not accessible by known 
renormalization group techniques. On the other hand, as we shall see, the equilibrium 
case A = 0, sometimes referred to as bounded Edward- Wilkinson, yields the simple 
rational result for the surface exponent 6ew = 3/4. In light of this result, it would 
be desirable to know if also the nonequilibrium surface exponents have simple rational 
values, an indication that the surface critical behavior is a trivial consequence of the 






(7) 
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known KPZ bulk properties. In particular, in Ref. [20] the value 6 = 7/6 has been 
conjectured for A < 0. Moreover, it is not yet fully clear how the contact points between 
interface and substrate are distributed and how finite-size effects influence the dynamics 
at the substrate. 

As a first step toward a better comprehension of the bounded KPZ class, we present 
here a detailed analysis of a simple solid on solid model generally belonging to the 
KPZ universality class, the so-called single step model (SSM). Once provided with an 
additional boundary restriction without additional short-range forces, it shows a simple 
wetting transition belonging to the bounded KPZ universality class [21]. In particular, 
the simplicity of the model allows for an analytical knowledge of it's wetting transition 
critical point = 0. In Sect. El the efficiency of its algorithm allows us to improve 
current numerical estimates of the exponent 9. 

In Sect. 0] we investigate the distribution of intervals between the contact points, 
identifying significantly different scaling behaviors depending on the sign of A. Finally, 
in Sect. EJ we study the scaling properties of finite systems, analyzing the time that 
is needed before the interface detaches globally for the first time. This kind of finite 
size scaling analysis is a common tool in numerical studies of nonequilibrium systems 
with an absorbing state, such as Directed Percolation (DP) [22], where the average 
absorbing time is known to scale according to the dynamical exponent, r ~ L z . 
However, supported by numerical results and analytical arguments, we argue that in 
the case A < this detachment time does not grow as a power law with the system 
size. Conclusions are drawn in Sect. El while in Appendix A we analytically compute 
the surface exponent when detailed balance is imposed in the SSM, thus proving the 
equilibrium result 9ew = 3/4. 

2. Measurement of the surface exponent in the single step model 

In what follows let us consider the so-called single step model (SSM) in one spatial 
dimension, which is probably the simplest and most compelling lattice model for KPZ- 
type interface growth. In this model the growing interface is represented by a set of 
integer heights hi residing at the sites % = 1 . . . L of a one-dimensional lattice of length 
L with periodic boundary conditions, obeying the restriction 

h i+1 -hi = ±l. (8) 

The model is controlled by a single parameter p e [0, 1] and evolves by random sequential 
updates as follows: Choosing a random site % the height hi is increased by two units with 
probability p provided that the restriction (jHJ is not violated. As shown in Fig. El this 
move can be thought of as depositing a diamond. Similarly, the height hi is decreased 
by two units with probability 1—p provided that the restriction (jHJ) is not violated. As 
usual each local update attempt corresponds to a time increment of dt = 1/L. 
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Figure 2. Single step model in 1+1 dimensions. Left panel: The simulation 
starts with a flat interface in the form of a horizontal sawtooth pattern. On 
selecting a site at a local minimum a diamond (rhombus) is deposited with 
probability p, flipping up the interface by two units. Similarly, selecting a site 
at a local maximum, a diamond evaporates with probability 1 — p, flipping 
the interface downward by two units. Right panel: Roughening interface after 
several deposition and evaporation events. 



For p = 1/2 the propagation velocity of the interface is zero, detailed balance 
holds in the stationary state and its dynamics is effectively described by an Edwards- 
Wilkinson equation [23]. For p ^ 1/2, however, the SSM exhibits KPZ growth with A 
being proportional to ~ —p. In this case the propagation velocity is nonzero and depends 
on the roughness and the average slope of the interface. For example, the initial velocity 
of a flat interface (a sawtooth pattern as shown in the left panel of Fig. |2J)) is given by 

v(0) = 2p-l. (9) 

As time proceeds, the interface roughens and the propagation velocity decreases until it 
saturates at a certain value -^(oo) which depends on the system size. For the KPZ class 
the so-called excess velocity v^it) — v L (oo) is known to decay with time as t -2 / 3 [24]. 
As a major breakthrough, Prahofer and Spohn recently succeeded in computing the 
rescaled height profile of a roughening KPZ analytically [25]. 

In contrast to many other KPZ growth models the SSM allows one to calculate 
the asymptotic growth velocity vl(oo) exactly. To this end one identifies segments of 
positive (negative) slope with particles (vacancies), mapping the single step model onto 
a partially asymmetric exclusion process (ASEP) [26] of diffusing particles with density 
1/2. Since the ASEP is known to evolve towards an uncorrelated product state with a 
current j = p/2 — 1/4, one can easily show that the propagation velocity of the interface 
in a finite system of size L » 1 tends to 

v L {co) = Km v L (t) = (p- (l + . (10) 

In order to study nonequilibrium wetting we need to introduce a hard-core substrate. 
As there is no parameter in the SSM to control the asymptotic growth velocity of the 
interface independently from the KPZ nonlinear term A, the substrate itself has to move 
so that we can approach the wetting transition by tuning its velocity v s . Obviously, this 
transition takes place even in finite systems and the critical point is exactly given by 
v s = v L (oo). 
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Figure 3. Measurement of the exponent 9 in the single step model. The left panel 
shows the density of contact points p s (t) as a function of time for p=l (A<0), p=l/2 
(A=0), and p=0 (A>0). The right panel shows the same data multiplied with the 
estimated power law t 6 . In addition we plotted the curve for p=l (A<1) multiplied by 
i 7 / 6 , demonstrating that the exponent 9 = 7/6 conjectured by Droz and Lipowski is 
incompatible with the numerical data. 



To determine the surface exponent 9 we simulated the SSM for p = 1/2 as well 
as for p = 0,1, where the non-linear effects of the KPZ term are most pronounced. 
In order to minimize finite-size effects we use a very large lattice size L = 10 6 . For 
p = 1/2 (corresponding to the equilibrium case A = in the KPZ equation) the interface 
velocity at the transition vanishes so that the substrate can be implemented by rejecting 
all updates that would lead to a negative height. The process is then simulated as 
usual and the density of sites at zero height is averaged over many independent runs. 
For p = 1 (corresponding to A < 0) we advance the substrate in temporal intervals 
At = 2(1 — 1/L) by one unit, flipping all sites below the new base line upward. Finally, 
for p = (corresponding to A > 0) the substrate moves backward by one unit every 
interval At, rejecting all updates that would render a height below the actual base line. 
Since the moving baseline effectively introduces a discrete time step in the algorithm, 
the density of contact points is measured immediately after (before) moving the wall 
when p = 1 [p = 0). 

The numerical results are shown in Fig. EJ Simulating up to 8 x 10 5 time steps 
and averaging over 160 independent runs one obtains convincingly straight lines in a 
double-logarithmic plot. Extrapolating local slopes we arrive at the estimates 

f 0.228(5) forp = 
9 = I 0.750(5) for p= 1/2 (11) 
[ 1.184(10) for p = 1 
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In order to illustrate the influence of scaling corrections, we multiplied each data set 
with t e in the right panel of Fig. |3] so that the curves become horizontal for large 
t. Our results have to be compared to the early numerical estimates 9\ = 1.1(1) and 
9 2 = 0.215(15) for the MN1 and MN2 Langevin equations [18]. A more recent numerical 
study of the MN1 equation yielded 9 X = 1.21(3) [17]. 

The numerical result for p — 1 is particularly interesting as it seems to rule out an 
earlier conjecture by Droz and Lipowski [20]. Analyzing a synchronization transition 
between two coupled lattices of tent maps and making use of a hyperscaling relation 
between stationary and spreading exponents, they proposed that for A < the exponent 
9 should be given by the exact value 7/6 ~ 0.1666. With the accuracy of the present 
simulations, however, this value lies clearly outside the error margin. To illustrate this 
difference we plotted p s (t)t 7 ^ 6 in the right panel of Fig. EJ As can be seen, the curve 
has a negative slope over at least two decades before finite-size effects set in at t ~ 10 5 . 
Our result suggests that in both cases of A 7^ the exponent 9 is probably not a simple 
rational number, rather it seems to be irrational. For A = (p = 1/2), however, our 
numerical estimate is in fair agreement with 9 = 3/4. In fact, as shown in Appendix 
A this value can be proven analytically provided that the two-point C(£) function of 
contact points decays algebraically as £~ 9z . 

3. Gap distribution of an infinite system 

Let us now investigate how the contact points of a critical interface with the substrate 
are distributed in an infinite system at a given time. Here an important quantity is the 
gap distribution G(£,t) of distances I between neighboring contact points measured at 
time t. As Fig. ((H) suggests, this distribution depends significantly on the sign of A. 

Our numerical results are shown in Fig. |3J For p = (A > 0) one observes an 
ordinary algebraic decay with a cutoff at the correlation length £ ~ t l l z . This type of 
decay is the same is in other non-equilibrium critical phenomena as, for example, in DP. 
For p > 1/2 (A < 0), however, the curves first decay algebraically, then reach a plateau, 
and finally cross over to an exponential cutoff. 

The qualitative difference between these two cases is due to the presence of two 
different length scales, namely, the bulk correlation length £ and the average gap size I. 
These two length scales grow with time as 



any point distribution on a one-dimensional line [27]. Depending on the values of the 
exponents z and 9, two different scenarios are encountered: 

• For A > (p < 1/2) we have 9 < 1/z and thus the average distance between 
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Figure 4. Panels A-C show the gap distribution function depending on the gap 
size measured in the SSM for p = 0, 1/2, and 1 at fixed times t = 16,32,64, . . .8192. 
Panel D shows a collapse of the data for A < according to the scaling form l|14fl . 
Contrarily, panels E and F demonstrate that this collapse does not work for A < 0. 



contact points £(t) is smaller than the correlation length £(t). In this case the gap 
distribution is a simple power law with an exponential cutoff at As in other 
critical phenomena this suggests the usual scaling form 

G(i,t) ~ t- {2 - ez) f y^y (14) 

In fact, plotting G(£,t)£ 2 ~ ez versus £/t 1 ^ z one obtains a convincing data collapse, 
as shown in Fig. HJD. 



• For A < (p > 1/2), where 9 > 1/z, the average distance between contact points 
is larger than the correlation length. In this case the gap distribution function is 
found to decay initially as £~ 0Z . Reaching the bulk correlation length I m £(t) it 
crosses over to an exponential function. However, in contrast to the previous case 
the associated length scale of this exponential function is £ ~ t e , hence in a double- 
logarithmic graph the curves exhibit a plateau at level oc t~ e extending from £(£) 
to £(t), followed by an exponential cutoff. 

In the following Section we will use these postulated scaling forms for the gap 
distribution to derive the first detachment time in a finite system, showing that the 
two cases lead to very different types of finite-size scaling. We note that the power-law 
scaling of the gap distribution is a postulate, supported by numerical results shown in 
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Fig. |U and the generally accepted power-law scaling of p s (t) = l/£(t) ~ t~ e . However, 
to our knowledge a rigorous proof of power-law scaling in the gap distribution is still 
missing. 

4. First detachment in a finite system 

We finally turn our attention to finite size scaling analysis, which a common tool to 
analyze properties of systems at criticality. We recall that the correct order parameter 
describing the wetting transition is the density p s (t) of sites pinned to the substrate, as 
defined in Eq. (JIJ). Comparison with DP naively leads one to draw an analogy between 
the pinned sites and the infected sites of systems with an absorbing state (i.e. a state 
in which the system dynamics get trapped forever). However, in the present case the 
depinned phase (i.e. the phase in which the substrate is completely covered) is not 
absorbing in the DP sense. In fact, since the wetting transition takes place at zero 
relative velocity, fluctuations can easily bring the interface back to the substrate, thus 
leaving the depinned state. Hence, instead of studying the average survival time r abs 
(i.e., the average time at which the system becomes trapped in an absorbing state), here 
we are forced to consider merely the first depinning time 77, i.e., the average time at 
which an initially exposed substrate gets completely covered for the first time. 

While in systems with an absorbing phase one expects the average survival time in a 
system of size L to scale as T a b s ~ L z , where z is the dynamical exponent, our numerical 
results will show an apparent power law behavior 77- ~ L 7 , where the exponent 7 depends 
on the sign of A and is generally different from the KPZ dynamical exponent. However, 
a more detailed analysis of the first depinning time probability distribution Pl(£) will 
show that at least for A < the scaling form 



is not satisfied, thus indicating that the observed power-law behavior does not hold 
asymptotically. 

4-1. Numerical results 

In order to measure the first depinning time 77, we used the single step model with a 
moving hard-core substrate as described above. To test finite size scaling we simulated 
various system sizes ranging from L = 64 to L = 16384. 

For each realization we determined the time t at which the interface was completely 
depinned for the first time. The resulting histograms of the probability distribution Pz,(i) 
of first depinning times are displayed in Fig. Moreover, we calculated the average 
depinning time for different system sizes and the relative width of the characteristic 




(15) 




Figure 5. Panels A-C show the distributions of the first depinning times for for 
p=0, 0.5 and 1. The system size increases from left to right starting with L=64 and 
going up to L=1024, 2048 and 16384 respectively. In panels D-F we test a collapse of 
the the above distributions according to Eq. lfT5|l . Only for p — the collapse works 
convincingly. 



distributions 

At* 

Orel = (16) 

where Ary is the standard deviation of first depinning times. The dependence of both 
quantities on L is shown in Fig. El Our numerical results can be summarized as follows: 

• In the case of A > we simulated the SSM with p = (i.e., no absorption just des- 
orption), averaging over 1.56 x 10 6 realizations for each system sizes up to L = 1024. 
We find that there is a strong indication of an ordinary power law behavior. The 
average depinning time scales with L 7 , with 7 = 1.41(2), while the probability dis- 
tributions Pl{£) have a constant relative width and, as it is shown in Fig. El (panel 
D), they can be cleanly collapsed according to Eq. (|15j) . 



• To study the case A = and A < we simulated the SSM with p = 0.5 and p = 1 for 
system sizes ranging from L = 64 to L = 2048 (p = 0.5) or to L = 16384 (p = 1), 
averaging over 8 x 10 5 realizations. In these cases too, the average depinning time 
Tf shows an apparent power law scaling. However, the relative width a re i of the 
probability distributions shrinks with increasing system size so that the probability 
distribution curves cannot be convincingly collapsed, making it obvious that we 
have no ordinary scaling behavior (see Figs. El and EJ). As it will be shown in 
the next section, this rather unconventional behavior can be traced back to the 
particularities of the gap distribution. 
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Figure 6. Left: Mean depinning time r/ plotted against the system size L. Right: 
Relative width cr rc i of the distributions as a function of L. As can be seen, the relative 
width stays constant for p = 0, in contrast to the cases p = 0.5 and p = 1. The solid 
lines represent the theoretical predictions for the cases p = 0.5 (blue) and p — 1 (green) 
according to Eq. (J22J)- 

4-2. Analytical arguments for the case A < 

In the following we show how the unconventional behavior observed for A < can be 
explained analytically in terms of the gap distribution function. The calculation is based 
on the observation that for A < the moment of first detachment in a finite system 
typically takes place long before the correlation length £ (t) reaches the system size. This 
means that the interface may be thought of as being divided into L/C,(t) uncorrelated 
segments, each of them touching the substrate with a small probability (cf. Fig. Ej). 
The moment of first detachment is then interpreted as the first fluctuation in which all 
these segments happen to detach independently. This assumption justifies the following 
approximations: 

(a) The distribution of first detachment times in a finite sample of size L with periodic 
boundary conditions approximately coincides with the corresponding distribution 
of first detachment for an arbitrarily chosen segment of size L in an infinite system. 

(b) The distributions of contact points at different times are assumed to be 
uncorrelated. 

These two assumptions allow us to use the gap distribution investigated in Sect. El and 
to compute the distribution of first detachment times analytically. While assumption 
(a) is reasonable in systems much larger than the correlation length, assumption (b) is 
crucial as it neglects all temporal correlations of contact points. 
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Figure 7. Schematic sketch of a critical interface in a finite system for A < at a fixed 
time t, divided into mutually uncorrelated segments of size The contact points 

are marked as red dots. The moment of first detachment corresponds to a fluctuation 
where all segments independently detach from the surface for the first time. 



Obviously, the probability for a segment of size L in an infinite system to detach 
at time t is given by 

ndm f?di(e-L)G(e,t) 

LjL[t) ~ f™dt£G(£,t) ' 1 ' 

where G(£,t) is the gap distribution discussed in Sect. El Consequently, the probability 
that the segment never detached up to time t is given by the product of l—Q d L (t) dt 
over all time steps or - in a continuous formulation - by an exponentiated integral 



<&(*)= exp(-* f Q dt'Q d L {t'] 



where k is a constant to be fitted to the numerical data. The quantity of interest, 
namely, the probability of first detachment, is then given by 

P L (t) = - j t Ql(t) . (19) 

An examination of these expression reveals that the probability of first detachment is 
determined by the occasional occurrence of very large gaps. This means that Pi{t) 
is governed by the right edge of the plateau shown in Figs. I3J5-IIP while the initial 
power law decay of the gap distribution has virtually no influence. Therefore, we may 
approximate the gap distribution function by 

G(t,t)~At- 9 exp(-B^J (20) 

where A and B are non-universal constants. With this "exponential" approximation 
the above integrals can be solved, giving 

Q d L (t) = exp (-B-) (21) 



and 



P L (t) = k exp (-Bj e - k{B Q )1/6 r(-l/6, BL/t e )\ , (22) 
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Simulation p=0.5 L=64 
Simulation p=0.5 L=1024 
Theory p=0.5 L=64 
Theory p=0.5L=1024 





Simulation p-1 L-64 
Simulation p-1 L=16348 
Theory p- 1 L=64 
Theory p=l L=16348 




Figure 8. Overlay of analytical predictions for the first depinning time probability 
distribution for A < (p=0.5 and 1) and numerical results of our simulation on selected 
system sizes. 



where T(a,z) = f^°t a ~ 1 e t dt is the incomplete gamma function. Here the scaling- 
invariant combination L/t e appears in two places whereas the prefactor of the gamma 
function is proportional to L 1 ^ . Obviously it is this prefactor which breaks scaling 
invariance so that Pi{t) does not satisfy the scaling form (jl5j) . 

In Fig. El the average depinning time and the relative width computed numerically 
from Eq. (|22j) for p — 1, 1/2 are found to be in good agreement with numerical data, 
while Fig. |H] shows a direct comparison between the analytical and the numerically 
derived probability distribution functions for the first depinning times. The parameters 
B and k have been determined by fitting to the numerical data and to grant 
normalization. Discrepancies between theory and simulations mainly occur for small 
system sizes and can be traced back to the above approximation of the gap distribution 
by a exponential function in Eq. ()18j) . As can be seen in Figs. the initial power- 

law part of the gap distribution function neglected in the theory is expected to play a 
more pronounced role for small system sizes. However, for large system sizes, where the 
"exponential" approximation used to derive Eq. (j2~Tf performs better, the agreement 
steadily improves. 



5. Conclusions 

Nonequilibrium wetting can be modeled by a KPZ growth process combined with a 
hard-core substrate that gives rise to an additional potential term in the KPZ equation. 
Adding such a term breaks up-down symmetry and leads to non-trivial surface features, 
which depend on the sign of the parameter A of the KPZ nonlinear term. For numerical 
studies of nonequilibrium wetting transitions we advocate the use of the single-step 
model because of its simplicity and efficency 

With large-scale simulations presented in this paper we have improved numerical 
estimates of the surface exponent 6. These estimates strongly suggest that 9 takes non- 
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trivial irrational values for A ^ 0. In particular we can rule out the simple rational 
value 9 = 7/6 for A < suggested by Droz and Lipowski. The failure of their argument 
can probably be traced back to the absence of an absorbing state in the bounded KPZ 
transition, questioning the feasibility of spreading analysis and hyperscaling relations 
used in absorbing phase transitions in the present case. However, in the equilibrium 
case (A = 0) our result confirm the simple rational value 9 = 3/4, in agreement with 
analytical arguments. 

In order to investigate finite-size scaling of wetting processes described by the 
bounded KPZ equation we introduced the average first depinning time as opposed to 
the average survival time in absorbing phase transitions. In one spatial dimension this 
naturally leads one to consider the distribution of gaps between pinned sites. We find 
that two different length scales are shaping the surface dynamics, namely, the spatial 
correlation length £(£) ~ t l / z and the average gap length l(t) ~ t e . If this latter length 
is larger then as in the case A < 0, gaps between contact points effectively behave 
as uncorrelated, thus breaking scale invariance for the first depinning times. This result 
comes as a surprise if one naively considers the analogy with absorbing phase transitions, 
where the scaling between temporal and spatial quantities is simply dictated by the 
dynamical exponent z. 

The case of A > on the other hand, shows a clean finite size scaling behavior 
of the first depinning times. Since in this latter case the spatial correlation length is 
larger than the average gap size, the analytical arguments presented in Sect. 14.21 no 
longer apply so that we are not able to predict the shape of the first depinning time 
probability distribution Pl(£). However, we note that our numerical result 7 = 1.41(2) 
is less then 10% off the known value for the KPZ dynamical exponent z = 1.5. Since 
it has been proven in [18] that spatial and temporal correlations are not asymptotically 
altered by the additional boundary condition, we expect 7 to converge gradually versus 
z in the limit L — > 00. However, for large but finite systems within the numerically 
accessible range the estimates may deviate. 
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Appendix A. Derivation of the exponent 9 in the equilibrium case 

In this Appendix we derive the exponent 9 for the equilibrium case p = 1/2 in the 
SSM analytically. In this case all transition rates between mutually reachable interface 
configurations are equal. The system thus approaches a stationary equilibrium state 
obeying detailed balance in which all interface configurations compatible with the 
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restrictions \h{ — h i+ i\ = 1 and hi > occur with the same probability. Although in 
an infinite system this state is never reached, a well-defined quantity is the conditional 
correlation function 

C (£) = ( ^;°/ fe 'r ,o) (A.i) 

where (■) denotes the spatial average over the index i. This correlation function may be 
interpreted as the conditional probability to find site i + £ at zero height given that site 
i is a contact point, too. As the interface roughens, numerator and denominator tend 
to zero but their quotient tends to a well-defined non-zero value. 

The present proof relies on the assumption that this two-point function - like other 
two-point functions in scale-free situations - decays according to the power law 

c{£) ~ £~ 6z , (A.2) 

where z = 2 is the dynamical exponent. As will be shown, the exponent 9 can be 
computed by deriving an exact expression for c(£). 

Adopting a transfer matrix approach used in Ref. [28, 29], the correlation 
function (|A.ljl may be expressed as 

(l\T e \l) 

C{£) = , (A.3) 

where T is an infinite-dimensional transfer matrix with matrix elements 

T *Ho otLSr 1 fc,V 6 {0,1,2,...} (A.4, 

while A = — 2 its dominating eigenvalue and (1| = |l)t = (1, 0, 0, . . .) denotes the vector 
representing height zero. Because of the special structure of T, where only the secondary 
diagonals are occupied, the evaluation of Eq. (jA.3|) reduces to counting all possible paths 
of the interface from site % to site % + £. The number of such paths is given by Catalan 
numbers, leading to the exact solution 

if £ even 



C(£) = \ (WW2TW (A.5) 
^ if £ odd. v ' 



Inserting Stirling's formula n\ ~ \/2ixn n+l l 2 e n it straight forward to show that this 
correlation function decays asymptotically as 

c{£) ~ J-e£- 3/2 (A.6) 

V 71 



Together with the known exponent z = 2 we thus arrive at 9 = 3/4. 
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